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ABSTRACT 


The cross section for photoneutron production in Be” at energies of 
1.692, 1.720, and 1.78 Mev has been investigated using radioisotope 
gamma-ray sources. Neutron counting was by manganous sulfate bath 
and a long counter. Gamma-ray counting was by a scintillation spec- 
trometer consisting of a Nal(T1) scintillation crystal and a 256-channel 


pulse+height analyzer. The following results were obtained: 


01 692 * 1.35 + 0.08 mb 


Oi. 70 
01. 692 


= 0.763 + 0.022 


ees 
OT. 692 


= 0.617 + 0.013 


A theoretical treatment of the P~>S transition near threshold using 


square=well potentials is also presented. 
I. INTRODUCTION 


The photoneutron cross section in beryllium is of great interest for 
several reasons: (1) beryllium has the lowest threshold for this reaction 
of any element; (2) beryllium is used in many photoneutron sources where 


a supply of monoenergetic neutrons is desired; (3) beryllium is often used 
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as a reflector in nuclear reactors; (4) the photodisintegration cross sec-~- 
tion may be used as a test for models of the beryllium nucleus. The pho- 
toneutron reaction was discovered in 1934 by Chadwick and Goldhaber 
(Ch34), who observed the reaction in deuterium, and in the same year 
Szilard and Chalmers (Sz34) observed the reaction in beryllium. Early 
measurements of the cross section for photoneutron production in beryl- 
lium were attempted by several experimenters, with widely varying re- 
sults. The cross section was calculated theoretically by Guth and Mullin 
(Gu49) in 1949, and their results have been fairly well accepted until very 
recently. The most recent theory is essentially similar to that of Guth 
and Mullin, but substitutes a potential of the Saxon type for the square well 
Originally used, resulting in a considerably better fit with experimentally 
determined points. 

This paper will present a brief discussion of photoneutron reactions 
in general, some discussion of the theory, a discription of several new 
experiments measuring the cross section at different energies, and con- 


clusions drawn from these experiments. 


Il. GENERAL 


Neutrons were first discovered in 1930, having been produced by the 
(a,n) reaction. The (y, n) reaction was then discovered in 1934, when 
neutrons were observed to result from the irradiation of beryllium and 
deuterium by high-energy gamma rays. Further investigation showed 
that these photoneutrons were emitted with a well-defined energy spectrum, 
in contrast with those resulting from alpha-particle bombardment, which 
have a widely distributed energy spectrum. In many applications, this 
discrete energy characteristic is very useful. Unfortunately, however, 
there is an attendant disadvantage. The cross section for the production 
of photoneutrons is very small, so that in order to produce a reasonable 
neutron flux a gamma-ray source of very high intensity must be employed. 
Unlike the alpha particles, which are stopped within the source, the gamma- 
ray intensity is attenuated only negligibly, thus producing high gamma-ray 


backgraunds and making handling conditions difficult. 
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(y,n) reactions in general involve the absorption of a gamma photon 
into the nucleus, forming an excited nucleus which then decays with the 


emission of a neutron. Thus 


EW = aA — (binding energy). 


From this it becomes apparent that if the gamma-ray energy is well- 
defined, the neutron energy will also be well-defined. The best source 
of monoenergetic gamma rays is the radiation from radioisotopes. Un- 
fortunately, the useful gamma rays (that is, those having a long half-life 
and high yield) are limited to energies of less than 3 Mev. Although the 
(y,n) reaction occurs in many eenents, only deuterium and beryllium 
have thresholds below 3 Mev for this reaction. By the use of various 
radioisotopes these two elements have been employed in sources for the 
production of monoenergetic neutrons with energies from 25 to 1000 kev. 
Photoneutrons can also be produced by the use of bremsstrahlung from 
targets bombarded by high-energy electrons from an accelerator, but 
these neutrons are not monoenergetic due to the continuous energy spec- 


trum of the bremsstrahlung. 


Ill. THE (y,n) REACTION IN BERYLLIUM 


Photoneutrons are produced from beryllium by the following reaction 
for gamma-ray energies less than about 2.5 Mev: 


7 


hv + Be Be n—1.666 Mev 


In this reaction the ee is left in its ground state and decays with a half- 
life of ~10~ > second as follows: 


Be > 2 Het + 0.09 Mev 


It would appear that the reaction 


hy + Be +2Het + n+ 1.576 Mev 


1s favored energetically, since it has a lower threshold. In 1952 Alburger 


(A152) looked for this reaction, using 1.631-Mev gamma rays from New 
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and observed no neutrons which could be attributed to this reaction. On 


the basis of this experiment he concluded that the cross section for the 


reaction could not be greater than 10 


=3 0) ed 
Cm. 


At gamma-ray energies higher than about 2.5 Mev many other re- 


actions are possible. 


cent paper (Pa60) and will not be repeated here. 


These have been summarized by Paolini ina re- 


The total cross section for the Be’ (y, nha: reaction has been meas- 


ured experimentally at several energies using radioisotope sources. 


Some of the more important of these are listed in Table I. 


Gamma 
source 


op l24 
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TAB ES! 


Measurements of Photoneutron Cross Sections for 


Beryllium Near Threshold 


Gamma energy 
(Mev) 


eo 7.2 


AA’) 


254,05 
Coo! 


2.754 


Unpublished. 


Cross section 


(mb) 
1.0mb 
Weg 
1.262 + 0.069 
[35 = 0706 
1.03 + 0.07 
0.83 + 0.06 
0.654 +4 0.031 
O55 
0.39 
0.50 
O23 77 + 0..02 
0.70 
0.674 





Reference 


Note 1 
Note 2 
(Gi59) 
This paper 
This paper 
This paper 
(Gi59) 
Note 3 

(Ha 53) 
(A150) + (Bi50) 
(Ed56) 
(Ru48) 

(Sn 50) 


Found by A. H. Snell about 1946. 


Obtained from (Ru48) by application of branching ratio. 


Unpublished. 


Found by E. Segré in 1949. 
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In 1956 Connors and Miller (Co56) explored the cross section near 
threshold using bremsstrahlung from an accelerator and reported a peak 
at 1.70 + 0.05 Mev, with the cross section > 200 times the average cross 
section in this energy region, and with a peak width < 50 kev. A more re- 
cent experiment by Walter, Shea, and Miller (Wa60) fails to show sucha 
peak. In this experiment they also used bremsstrahlung from an electron 
accelerator, using a thin gold target, and varied the maximum electron 
energy in 5-kev steps. Their results agree with the cross sectionas 
measured with radioisotope sources. They analyzed the bremsstrahlung 
spectrum using the Sauter-Fano theory, finding that this gave a better fit 
than did the earlier Bethe-Heitler theory. They attribute the peak reported 
in the earlier experiment to an incorrect analysis of the bremsstrahlung 
yield from the thick gold target used at that time. 

A recent experiment at M.I.T. by Paolini (Pa60) determined the over- 
all photoneutron spectrum from about 3 to 18 Mev. This is shown in Fig. 1. 
In this figure the peak just above threshold as determined with monoener- 


getic sources has been added to Paolini's results. 


IV. Tere ORY 


In 1949 Guth and Mullin (Gu49) calculated a theoretical cross section 
for photoneutron production in Ber. based on a two-body model consisting 
of a single neutron loosely bound toa Bee core, and using square-well 
potentials. In their model the ground state is takenasa Peay 2 state. Elec- 
tric dipole transitions then lead to an Ss) /2 state near threshold and to Saye 
and /2 states at somewhat higher energies. The cross section fora 
magnetic dipole transition, which would lead toa emt / 2 state near thresh- 
old, was calculated to be very small. Their expression for the cross sec- 
tion for the P — S transition contains four adjustable parameters namely, 
the radius and depth of the initial and final wells. Making the assumption 
that the initial and final wells have the same radius of action and choosing 
the initial well depth to give the correct binding energy, reduces the num- 
ber of adjustable parameters to two. For a given radius the depth of the 


final well is then chosen to result in a fit to the experimental cross section 
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near threshold. In their calculations they used a radius of 5.0 f (1 fermi 

1@ 42 cm) and a final well depth of 3.0 Mev. The threshold (binding 
energy) at that time was thought to be 1.63 Mev. The total cross section 
calculated by Guth and Mullin is shown in Fig. 2. 

An early theoretical treatment by Bergmann (Be51l) was similar to 
that of Guth and Mullin, but the numerical results were different. Some 
of the data included in his paper indicated the effect of varying the pa- 
Parie ters. 

It appeared to us that the use of more recent experimental data, to- 
gether with a judicious choice of parameters, might produce a better fit 
to the experimental cross section. On this basis, new square-well cal- 
culations following the treatment of Guth and Mullin were made. The de- 
tailed development of this theory is included in Appendix A of this report. 
Calculations were made using several values for the radius of the initial 
and final wells, and the initial well depth in each case was chosen to give 
the correct value of 1.666 Mev for the binding energy of the last neutron. 
The depth of the final well in each case was chosen for two different con- 
ditions: (1) correct energy dependence of the cross section, based on the 
experimental ratio of the cross section at 1.78 Mev to that at 1.692 Mev; 
and (2) correct value of the cross section at 1.692 Mev as experimentally 
determined. The results of these calculations are shown in Figs. 3 
through 5, which also illustrate the effect of changing the various pa- 
rameters. Itis felt that, because of the magnitude of the parameters 
used to fit the experimental cross section, the square-well assumption is 
not entirely valid. The qualitative fit, however, is quite good, especially 
considering the simplicity of the model. It is interesting to note that the 
sharp peak just above threshold demands that the final state be a lightly 
bound state (~50 kev binding). 

Goldman, Francis, and Guth (Go60) have recently calculated the 
theoretical cross section near threshold for the P —S transition using a 


single-particle shell model of the Saxon type: 
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r —— THEORETICAL CURVE 
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X SNELL,et al. 
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Ul oe Be’(y,n) cross section (from Gu48). 
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EXPERIMENTAL: 


@ JOHN & PROSSER 
© GIBBONS ef al. 


THEGRETICAL.  -2o7e 


( USING SQUARE WELLS) Vo, Yo= 9.06 Mev. 
REDUCED X0.713 

4 NORMALIZED. Vig 

, =2.59 
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0.2 

O 
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Jo oe Bec) cross section curves illustrating the effect of 
varying rg- Vo 2 is chosen to fit the ratio of the cross sections 
at 1.692 and 1. 78’Mev. Note that the absolute magnitudes have been 
multiplied by reduction factors. 
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-- pe sof | Vorve FOR CORRECT 
eee o= 


ABSOLUTE VALUE OF 
lo= |.35f O| 692 


EXPERIMENTAL. 
@ JOHN & PROSSER 


O GIBBONS ef al. 





1.66 1.70 .75 1.80 1.85 90 


MUL-10016 
Ey ( Mev ) 


Fig. 4. For three vajnes of Too Vo We has been chosen to fit the 
absolute value of the Be’(y,n) cross OMgon at 1.692 Mev. 
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— THEORETICAL, USING 
SQUARE WELL WITH 
lo =|.87f & Vo,1/2 
CHOSEN FOR CORRECT 91.692 


EXPERIMENTAL: 


@ JOHN & PROSSER 
O GIBBONS, et al. 


Ops (mb) 





1.66 .7O Neal a, 1.80 1.85 L-9@ 


MUL-10027 
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Figs 5. The Bene) cross section curve with ro and V 1/2 chosen 
to fit the absolute magnitude at 1.692 Mev and the ratio of Pad cross 
sections at 1.692 and 1.78 Mev. 
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—Vo 


r—-R 
a 


AC 
1 + exp 





with Vg determined from the binding energy of the last neutron, and with 
the final well based on scattering data. Using R = 3.0f for both the in- 

itial and final states, witha = 0.6f for the initial state anda = 1.2f for 
the final state, they have obtained a cross section curve which appears to 


fit the experimental data in both energy dependence and absolute magnitude. 


V. EXPERIMENTAL PROCEDURE 


A. Sa Gamma-Ray Source 


7 for the 


It was decided that the total photoneutron cross section of Be 
1.692-Mev gamma rays from yee would be measured. This was the 
closest gamma-eray energy to the threshold which was readily available in 
sufficient intensity and with a long enough half-life to permit convenient 
measurement. Also, this energy appears to be quite near the top of the 
peak in the cross section near threshold. Furthermore, So sources 
had been used in previous experiments (see Table I) which had resulted 
in cross section Measurements differing by almost a factor of 2, and it 
was desired to clear up this long-standing discrepancy. 

A source consisting of a hollow sphere of 1/16 in. aluminum packed 
with granulated antimony to a density of about 2.7 esol and having a 
1.36-cm outside diameter, was irradiated in the Livermore Pool Type 
Reactor (LPTR) to a total antimony activity of 225 millicuries and allowed 
to decay for a period of several months to eliminate any short-lived activity. 
This gamma-ray source was then held within a hollow beryllium sphere 
made up of nesting hemispheres, the whole assembly of beryllium having 
an outside diameter of 8.2cm anda thickness of 0.645cm. A sketch of 
the assembly is shown in Fig. 6. The particular source geometry used 
was chosen to: (1) give spherical symmetry; (2) minimize the correction 
necessary for the finite size of the Syl source; and (3) reduce the ab- 
sorption of gamma rays and the moderation of neutrons in the beryllium 


shell to a reasonable amount. 
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BERYLLIUM 





8.2 CM 


Fig. 6. Beryllium photoneutron source. wyyf-8456 
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B. Neutron Counting 


In order to measure the neutron flux, a Jong counter was constructed 
following the pattern of the Harwell 4 model (A155). The effective face of 
the counter was determined by using Po-Be, Pu-Be, Po-lLi, Sb-Be, and 
mock fission sources at various distances, and plotting (count cme 
vs distance from the actual face of the counter. The use of the long 
counter for absolute neutron flux measurements proved unsatisfactory, 
however, because of the different geometries of the various sources and 
an uncertainty as to the variation in efficiency of the counter over sucha 
broad range of neutron energies. 

As a better means of counting the total neutron flux, a manganous 
sulfate bath was prepared (Ge59). This method of determining neutron 
flux has the advantage that essentially all the neutrons are moderated and 
captured by the manganese, so that the resulting activation is independent 
of the neutron energy spectrum. A tank 44 in. in diameter and 29 in. deep 
was filled with a nearly saturated solution of manganous sulfate, using 
175 gallons of distilled water in which was dissolved 575 pounds of MnsSO, . 
A Victoreen 1B85 Geiger-Muller tube was attached to a probe of aluminum 
tubing and waterproofed for dip counting. A eae” button source was fixed 
ina standard geometry for periodic checks of the G-M tube counting rate. 
An NBS Ra-Bée source was used to calibrate the MnSO, bath. It was sus- 
pended in the center of the tank ina waterproof holder made of aluminum 
and the bath was activated to saturation. Since the half-life of Mee is 
2.58 hours, any activation over a period of a day or longer essentially 
saturates the bath. The G-M tube was then placed in the center of the 
bath and a series of counts taken over a period of several half-lives of 
the Min? to ensure that the proper activity was being counted. From the 
results of this series of counts, the count rate at the time of removal of 
the source from the bath was computed. From this initial counting rate 
and the total flux of the NBS source, the total flux of any other neutron 
source can be determined by immersing it in the bath and activating the 


bath to saturation. The ratio of the fluxes is then the same as the ratio 


of the initial counting rates. The initial counting rate obtained with the 
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NBS source has proved to be reproducible to within 1% in subsequent 
measurements. This procedure was followed in determining the total 


neutron flux produced by the antimony-beryllium source assembly. 


C. Strength of Gamma-Ray Source 


In order to compute the photoneutron cross section, it was necessary 
that the total flux of gamma-ray photons with an energy of 1.692 Mev be 
determined. This was accomplished by the use of a thallium-activated 
sodium-iodide scintillation crystal and a 256-channel pulse-height an- 
alyzer. Computations of the gamma-ray intensity were based on un- 
published curves for intrinsic crystal efficiency and photofraction fur- 
nished by Harry West of this laboratory. It is estimated that the ac- 
curacy of this method is better than5%. The scintillation counting tech- 
nique is felt to be the best method of determining the strength of the 
gamma-ray source in this experiment, since it is independent of the de- 
cay scheme of the source, and provides an automatic correction for self- 
absorption and Compton scattering in the source. In using other methods 
the strength of a particular gamma line must be calculated from a know- 
ledge of the decay scheme involved, and in many instances the decay 
scheme is not known with great accuracy. With the scintillator and 
pulse-height analyzer combination, only those gamma rays that have a 
specific energy are counted, and further, any gamma rays absorbed in 
the source are not counted, just as they are not seen by the beryllium 
surrounding the gamma-eray source. In addition, any gamma rays which 
are Compton scattered are not counted under the photopeak. This is de+ 
sirable, in the case of antimony, since almost all (99.5%) of those photons 
which have undergone Compton scattering are below the threshold for 
photodisintegration. Thus the gamma-ray source looks the same to the 
scintillator as it does to the beryllium, and corrections for absorption 
and Compton scattering are applied automatically, giving a more ac- 
curate evaluation of the gamma-ray flux available to produce photo- 


neutrons. 
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D. Computation of Cross Section 


Once the gamma-ray flux (at the required energy) and the neutron 
flux are known, the total cross section for photoneutron production at 


the particular energy involved can be computed from 


N A 
o=> n ae 
y <0 a Be 
where 
Oo = cross section in cm 
No = neutron flux in neutrons/sec 
Ap. = atomic weight of Bev 

N = gamma-ray flux in photons/sec 
Ny = Avogadro's number 
Pp. = density of pen in g/cm? 


Ar = thickness of Be’ shell in cm. 


E. Comparison Technique 


It was desired to explore the photoneutron cross section near thresh- 
old in as much detail as possible using radioisotopes, but the gamma-ray 
sources available in this energy region are limited to a few which have 
short half-lives and/or low intensities. With a half-life less than several 
hours or a total neutron flux less than about oe n/sec, the use of the Mnso, 
bath becomes subject to large uncertainties. Therefore, a new experi- 
mental arrangement providing for the following considerations was de- 
sirable: 

1. The measurement should be independent of half-life. This could 
be accomplished by a simultaneous count of gamma rays and neutrons. 

2. It should be possible to count a weak neutron flux. This could be 
done with a long counter. 

3. The effect of room scattered neutrons should be minimized. This 
could be done by placing the source close to the long counter. 

4. Since there are about al gamma photons for each neutron prod- 


uced, and since the scintillation counter is more efficient than the long 
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counter, precautions must be taken to prevent excessive gamma count- 
ing rates, with attendant dead-time losses. This requirement could be 
satisfied by locating the scintillator at a distance from the source, and 
by using a collimator at the crystal as necessary. 

The arrangement finally devised is shownin Fig. 7. The simultane- 
ous count of neutrons and gammas results in an n/y ratio which is pro- 
portional to the cross section. When several gamma sources of differ- 
ent energies are used successively in the same geometry, a relative cross 
section curve is readily obtained. This can then be converted to an abso- 
lute cross section curve if the correct aboslute value is known at one of 
the energies. 

In computing the n/y ratios, corrections are made for the variation 
in intrinsic efficiency and photofraction of the Nal scintillator. No cor- 
rections were made for variations in efficiency of the long counter. This 
counter was of the Harwell 4 type which Allen found to be flat down to 
Sb-Be photoneutron energy (A155). Even for the typical 10-15% droop of 
the efficiency of a long counter, our efficiency uncertainty would not be 
greater than about 2%. The long counter was biased to eliminate y count- 
ing. 

This comparison technique was used to make cross-section measure- 
ments at 1.78 Mev with pe as the gamma-ray source, and at 1.72 Mev 
with Bee A sphere of solid aluminum the same size as the antimony 
source was activated in the LPTR and placed in the center of the beryl- 
lium shell. It was necessary to work rapidly, since the half-life of ee 
is only 2.3 minutes. This was accomplished by removing the source from 
the reactor by a pneumatic tube delivery system, so that the total elapsed 
time after removal from the reactor flux until the count of the Al-Be as- 
sembly was started was less than one half-life. A lead collimator having 
a 3/4-in. aperture was used at the scintillation crystal to cut down the 
counting rate for the gamma rays while still maintaining a good neutron 
counting rate. Several successive runs were made as the source decayed, 
and those runs having good counting rates without excessive dead-time 


corrections were used in computing the n/y ratio. The n/y ratio for 
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antimony in the same geometry was measured and the cross section at 


1.78 Mev computed by 
(n/Y) 78 


(ee? lee 


It should be noted that the ratio of the cross sections is considerably 
more accurate than the absolute value. 

For the measurement at 1.72 Mev a foil of high-purity lead was bom- 
barded in the 60-inch cyclotron at Berkeley to produce the reaction 
ppe?6 (d, 2n) ee The procedure was then the same as that used in 
the case of the aluminum source, except that the bismuth source was 
considerably weaker, so the collimator was not used at the scintillation 
crystal. The ratio n/y for the bismuth source was somewhat more dif- 
ficult to determine than for either the aluminum or the antimony, since 
some Bie was present from the Be a sae reaction. A rather 
complex gamma-ray spectrum resulted from this contamination, and it 
was necessary to separate out its effect on the basis of pulse height and 
half-life. 

There are several other radioisotopes with gamma-sray energies 
above threshold and short half-lives which could be used in this compar- 
ison technique. We know of no others, however, which have gamma-eray 
energies between threshold and the 1. 78 Mev from the Aneee Use of 


higher energies than this would probably demand an individual calibration 


of the long counter efficiency curve. 


VI. EXPERIMENTAL RESULTS 


In measuring the photoneutron cross section at 1.692 Mev, using 
gamma rays from S5 ae the following results were obtained: 
Strength of NBS Ra-Be source E-1277 as of 9 December 1954: 
Qo = (6.09 + 3%) X 10° n/sec. 
Strength as of 17 September 1959 considering growth of Po 
Q = Q, |1.12-0.12 exp ( -t/28)| 
(6.21 4.3%) X 10° n/sec. 


A110). 
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(dN/dt)) for NBS source in MnSO, bath = 1292 c/min 
(dN/dt)) for Sb-Be source assembly = 11.05 c/min 
N from Sb-Be = or KO. 2x 10° = 5.31x 10+ n/sec 
a for 1.692-Mev line of gp 2 as of 3 Sut ee 1959 = 25.602. a 
Me as of 17 September 1959 = 5.61 xX ae as (—0. 693 X 14/60. 9) 
= 4,78 10 photons/sec 


The following corrections and uncertainties apply: 





Correction to cross section Correction(%) Uncertainty(%) 
Gammaeray source strength --- Da 
Neutron counting statistics : --- 20 
Finite source size SU e <0") 
Attenuation of gamma ray in Be shell gene <0 
Effect of 2.09-Mev gamma ray —4.4 1530 
Neutron escape from MnsoO, bath —0.9 <0! 
Absolute strength of NBS source E-—1277 —— a. 0 
Total correction 3 il --- 
Overall uncertainty --- be 


The corrections for finite source size, attenuation of the gamma ray 
in the Be shell, effect of the 2.09-Mev gamma-ray, and neutron escape 
from the MnsSoO, bath are described in Appendix B. 

The cross section is now computed as follows: 

A 
og = NON, Pod ee O20) 
5.32 X 10* x 9. 02 X 0. 969 
4.78X 10° x 6.025 X 10°? x. 1.85 X 0. 645 
135 < 10e> teme 
1.35 + 07 0Semibre 


In determining the cross section at 1.78 Mev the comparison tech- 


O 


nique described previously was used, with Ace as the source of gamma 
rays, and the following results were obtained: 


n/y for Sb-Be = 0.1179 
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n/y for Al-Be = 0.0728 


Opens oe 1.35 = 0.83 mb. 

The ratio of the cross sections at 1. 78 Mev and 1.692 Mev is consider- 
ed to be much more accurate than the absolute value of the cross section, 
having an estimated uncertainty of + 2%. This gives an overall uncertainty 
at 1.78 Mev of 6.6%, so that 

0, 7g = 0-83 + 0.06 mb. vr 

In the case of the cross section at 1.72 Mev, the strength of the Bi 
gamma-ray source was rather low, so the collimator was removed from 
the sodium iodode crystal to obtain better gamma counting statistics. 

This change in geometry changed the n/y ratio for the Sb-Be assembly. 
The following results were obtained: 

n/y for Sb-Be = 0.0421 

n/y for Bi-Be = 0.0309. 

From these figures a preliminary cross section at 1.72 Mev was 
computed as 

ge i GOS, 


0) 9) = D0oarr * 1:39 = 0.99 eae 


It was necessary to make an additional correction to this prelimi- 
sanO5 , 
nary value, since there was an appreciable amount of Bi present in 


the source. pe has strong gamma rays at energies of 1.766 and 1.777 


Mev. These, together with the 1.72-Mev gamma ray from the Bio 
resulted in a broad photopeak made up of one peak at 1. 72 Mev and an- 
other, approximately one-fourth as intense, at about 1.77 Mev. Even 
if the ratio of the two peaks were 1:3 or 1:5 rather than the estimated 
ratio of 1:4, the resulting cross section would not change by more than 
about 1.3%. The assumption of a cross section of 0.85 mb as a reason- 
able value of 1.77 Mev leads to 

0. 80 72 + 0.2 (0. 85) = 0799, 
so that 

0, 7 * lal OS aig you 
Because of the complicated spectrum, an uncertainty of + 3% was 


estimated for the comparison procedure inthis case. This leads to 
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an Overall uncertainty of about 7% for the cross section at 1.72 Mev, 


with the result that 


01 7297 (i062) 0207) tabs 


In regard to the absolute value of the cross section, it should be 
mentioned that the antimony source was also counted by the Hazards 
Control Group using a Landswerk Condenser R-meter. The total 
strength of the source was measured as 47.5 mc as of 11] August 1959. 
The strength as of 17 September 1959 was then 31.2 mc. If we assume 
that the 1.692-Mev gamma ray occurs in 50.6% of the disintegrations 
(Gi59), the strength of this gamma ray was 5. 84 X 10° photons/sec. This 
iMmst be corrected for scattering within the source, since scattered 
gammas are counted by the R-meter but are below the photoneutron 
threshold. The correction is about 6.5%, so that 

Men] 5.46 lO enorone) coom 
Subsituting this value ane cross section formula we obtain 

Go = 1,176°2°0, 07 mob - 
As stated previously, it is felt that this method of measurement is not 
as reliable as that using the scintillation crystal. However, the difference 
in the values obtained by the two methods indicates the presence of some 
systematic error, possibly in the branching ratio. 

In an attempt to improve the accuracy of the original cross section 
measurement at 1.692 Mev, particularly the neutron counting statistics, 
a stronger antimony source has been prepared and the neutron flux 
measured. This source was also counted by the Hazards Control Group. 
Their measurement of the stronger source leads to the same cross sec- 
tion, within 1%, as did their measurement of the earlier source. The 
estimated uncertainty for the neutron counting statistics as shown in the 
table of corrections to the cross sections is based on the measurement 
of this stronger source. Even though the cross section based on scin- 
tillation spectrometry is considered more accurate than that based on 
R-meter measurements, the latter have proved consistent in successive 
measurements of an antimony source, within the known accuracy of the 


half-life. 
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VII. CONCLUSIONS 


The value obtained in the present experiment for the photoneutron 
cross section at 1.692 Mev is in reasonable agreement with the recent 
accurate determination by Gibbons et al. (Gi59). This agreement re- 
moves the long-standing large uncertainty in the value of the cross 
section at 1.692 Mev. 

The energy dependence of the cross section near threshold is now 
well defined by the combination of our results at 1.692, 1.720, and 1.78 
Mev with the data of Gibbons etal. (Gi59) and Walter et al. (Wa60). 

The use of a square-well potential in the theory of Guth and Mullin 
does not give a fit to the data unless the radius is made unrealistically 
small (1.87f). The new calculation by Goldman et al. (Go60) using a 
Saxon-Woods potential with reasonable parameters results ina good fit 


to the experimental data. 
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APPENDIX A 
THEORETICAL DEVELOPMENT USING 
SQUARE-WELL POTENTIALS 


This development follows that of Guth and Mullin (Gu48) but is 
presented in more detail. We considered it necessary to verify the 
formulae of Guth and Mullin because of the different numerical re- 
sults obtained by Bergmann (Be51). 


The general formula for the electric dipole photoneutron cross 


ss » M, (=) (1) 
ih 


section is 





in which 


M,(E) = Ze = | be dt. (2) 


u/m is the ratio of the reduced mass of the system to the mass of the 
residual nucleus, z = r cos @ is the electric dipole operator, and 

E = hv —¢, where e is the binding energy of the neutron to the residual 
nucleus. 

Using a ground P state, the selection rule Al = + 1 for electric 
dipole transitions permits P~>Sand P~>D transitions. Only the P>S 
transition will be treated here, since it is the one which is important 
near threshold. 


The expression for the electric dipole disintegration of Bee then 


a Sry € 2) (ureya 


where the summation extends over all non-vanishing matrix elements, 


becomes 





Z 
T 





2 (3) 


Assuming the spin interaction to be small and of a central (non-tensor) 
type, the eigenstates f,i are described by the quantum numbers E 


(energy), j, De andl. The eigenfunctions are 
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BT ae 


l | l 1/ 
w(l, jzelts, ee 17 (l1+=+m.) x 1 (6,9¢)a 





l 2 
Be Yim, + 5 109) Pn (4a) 
R, .(r) 
ae alg = I,J l a We 
~a¢iemy)!/2 y 1 (6,¢)8 (4b) 
Z j 1 aA PIP, 


The matrix element We b. dr of Eq. (3) is different from zero only if 
Al = +1, Aj = 0, + 1, and Am. = 0. an and Bo a: spin operators with 
the property that a es 8 B.S 1 and me Bo — a, = OF 

Fora aye initial state, ay can be + 3/2 eye Be i with equal prob- 
ability. Since the final By 2 state can only have she = + 1/2, a factor of 
1/2 in the final cross section results from the selection rule Am. = 0. 
Furthermore, Eq. (4b) would apply toa Pi /2 state, so that only Eq. (4a) 


need be considered in the case ofa P3/2 initial state. 


Considering m. = +1/2, the initial eigenfunction is 
EE TE 
7 ae 3/2 1/2 
For m, = —1/2, the initial eigenfunction is 
R (r) , ‘) 
Wy (1, By —1/2) = tie pay, _ (9: > ) a. 45 (21/2 Yy lO PP P| > (55) 


For m, = + 1/2, the final eigenfunction is 


i 
¥(, 1/2, #1/2) = Ro 1 j2(t) Yo 9 (09) | (5c) 
Por a, = —1/2, the final Bi NS As 
¥(0, 1/2, 1/2) = Ro 1 pple) ¥o, (6%) B| (5d) 


The expression (ue Z p. d7 can now be written as 


1/2 ‘5 * 7 ity. Ss | 
(2/3) (WV \Ro,1/2 ¥0,0 r cos @g (83/2 *1,0) r~ sin gdr d@d@¢ »& 
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Separating variables gives 


Ie) 2 * 3 %% 
(2/3) {Ro 1/2 So 3/2 r dr ({ Yo 0 wo 9 Cos g singd @dg¢. 


The normalized spherical harmonics are 
me Ne ; 31/2 
. (=) tip 


which results in 


2 
(2) * 3 Z 
a = One 2 Seay ae dr cos gsingd@gd®¢. 


Integrating the angle part gives 


1/2 3 
a) ("RG 0, 1/2 ee 


Considering the sales of 1/2 mentioned previously, the expression 


» (ope YW. dr 


2 
a je) R ic = 1 IR 
. Bee 








now becomes 





2 9 
where 
R = " p* R 3 d 
s =| Oo) onlay eae 
0 


This makes the expression for the cross section 








5 ere e* alee a 
— aa C ps 
By using hv = Mw andh = ae this becomes 
ofr 725) We : ps! 


The only remaining problem is the evaluation of Bg Using the 
spherical qquare-well potential, for the initial state 


VS ewe = 


Vas= 0 a 


and the radial wave functions are 


R Aid, (Br) ees Py 


A, /sin or cos 8) 
(Br) 


ey 2 
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7 l1+ar ek ee Eee 
are ay2 a. y : 
fe 
where ; 1/2 
a = [2 
2 
e 
2 
and sect: 
_ Fe 3/28 e) | 
| ; 8M 
with yp = reduced mass of neutron = 5 
Matching the logarithmic derivative = = atr = » gives 


— > sin er . - 
a 


¢ 
Be,rcetepr, = 1+ (1 + age (8) 


io] 


Normalizing by using { eho anal node sl gives 
0 b) 
2 ‘ \4 
Ay z 11 ‘He war, ) (8 } 


Zz ee Br, 
+(+ar,) as —]; ———>/=1 


2 
(Br) 


For the final ye state, assuming the same radius of action for the 


initial and final states 
= — < 
Mi MG 1/2 a 
Viewers. 0 aed 85 


The radial wave functions for this state are then 


= } < 
Ro 1/2 ae (yr) ae, 
BON sin yr 
0 sis 
_ Aten lite < = ) + 5, Pe 
Ro 1/2 > ?o SS J ore : 
. 2 
where y= |= (V, 1/2 1p Sy) 
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b/ 2 
and k (25-5) 
N 2 
Hf 
It will be noted that the wave function for r> r is that for a free 
particle with phase shift 6 
Matching the logarithmic derivative at r = Re leads to 


rev sin \/ ag iB, sin Ds 


A, cos a A 5, cos 65 


The normalization in this case is not so straightforward as for the 


initial state. It is accomplished as follows: 


{ (Ry ip)” r“dr = 1 (8) 
0 
vo © sin{k(r—r ) + 6 a 
C ° 0 eS 
Re 7 Jo we r- dr + BS { ao r dr=l. 
0 ee 0 & a 


The first integral is readily integrated: 


bay, Zz 
2 ° | sin yr : Z . 
A { Vr Ip. Gig I lyr, — cos Yr, sin Wil : 
0 0 2Y 
To integrate the second integral, assume for the moment a small energy 
spread such that dk = (ky —k 5) << 1 and write the integral 
= + § = + 6 
sin ea a r : | sin|k, (r ae ) alt? 


a rgiietali ge 
kor , 


For the upper Phas substitute g and rewrite as 


BS lim & ; 
| : | ~~ | 
RGRS 6 acore ae eed a eR Mies ee is 


Now substitute x = r— ae giving 


Z 
BY lim > y ri ve Neb 
ae ; ee iy sin’k x 65) sin (k, x + 6) : 
0 
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This can now be integrated to give 
zZ pee 

Bo lim {sin (kK, —k,) g 

ee es 

ky 5 goo Z ky k,, | 

pin cs + ks) g + 26, 
i 
Z (k, + KS) 

Now, since (Kk, —k,) << 1, all terms resulting from the left side of Eq. (8) 


are negligible except 


Sees: - Ate eevee 
k = ~~ 2(k, —k_) 
Sr aia (5) ge 
This can be rewritten in terms of the Dirac delta function as follows 


(see L. I. Schiff, Quantum Mechanics, McGraw-Hill, 1959, p. 51): 


2 
a a 
k ok Zz ] p} 





tai 
For normalization per unit dk interval then 
BC Eis 
me Z 
or By = K(2/m)i/2 
For use in Eq. (7), however, | a é must be normalized per unit 


energy interval. This is accomplished by multiplying B, by 


MZ 1/2 
(3) +3) 


dk\1/2 (4 22 
y T 


VEN e 8 Oy 


so that 


In order to determine BS it is now only necessary to subsitute the - 
normalized wave functions into Eq. (6) and integrate. This is a lengthy 


but straightforward integration, resulting in 
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by 2 
B, a? [(2/m) (uk/)] (2 + 2 (B°/a7) (+ ary) —4e ly? = BY) 


R = 
Psy (k2 + ye mot ve, 72 | (2 = Bz) 
ar (2 tar.) +k“ ro 1h ot | ty? —p*) /a* | CUR eG ne) 
+ 3.3 | 
(GL ae Vy ee 
3+ ar, + fees} (l+ar_) 
Tr 5 om) 2 Vie GOL yz, : 
(a +k’) 4 


APPENDIX B 


CORRECTIONS TO THE CROSS SECTION 


It is necessary to make several corrections to the cross section 


computed from the basic formula 


_ ae Be 


o = TN Neo oar 
aS 0 PBe ™ 


These will be described in the following paragraphs. 


1. Finite source size. We follow the derivation given by Gibbons 


etal. (Gi59) , 


a 
Pay ; 
red oo 
_ 2 ee 
ro 9 = R + r°—2Rrcos@ 
\ / R | 
\ bo . : 
x 
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Assuming a homogeneous spherical source of strength Q and radius 





a, the strength per unit volume is sO. 
4a 
If Nt = number of eocrey ama in the Be shell, then 
Zz 
dn = Ntgo ee = dV 
2 3 
p 4tmTa 


gives the number of neutrons per second due to source element dV. The 
total neutron flux produced is nf where n, = NtoQ is the result if the 
source were concentrated ina point at the center of the shell, and fis 


the finite source size correction givenby: 


-8 (fe. 38 SR 
 “Vieeee 
p 


“sin gdr dod 
eS fia — 2Rr cos @ 











This integrates to 


po GEM RS set 5. (Bate 
2 R-a 
Za 


Since a<<R, the logarithm may be expanded: 





= x— 1 
aa 2st +3 " Ley areal |. Bs 0) 
Letting x = pit ; = = os 

ae Ok pte | 


Subsituting in the equation for f: 


4 6 
3 7 2 2 a la 
f=~-——/f/-at=— —te=- + 4 
Bae Ne 15 2 5 pt | 
ee | 
5R 


For extreme accuracy this should be multiplied by the average value of 
sec a, but this correction is negligible when a << R. 


2. Attenuation of gamma ray in Be shell. 
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liege 
1. 2 


3. Effect of the 2.09-Mev gamma ray. 

In the case of the: Sb-Be source assembly, the neutron flux measured 
in the manganous sulfate bath includes not only the neutrons produced by 
the 1.692-Mev gamma ray, but also those produced by the 2.09-Mev 
gamma ray. A correction for this effect must be made in computing the 
cross section at 1.692 Mev. The ratio of the 2.09-Mev gamma ray to 
the 1.69-Mev gamma ray was measured as 0.148 + 0.008 by Nal count- 
ing. A cross section of 0.4 mb was assumed as a reasonable value at 
2.09 Mev. The correction to be applied to the cross section is then 

0.148 X 0.4 = 0.06 mb. 
This was divided by the computed cross section to give a percentage 
correction. Our correction is the same as that applied by Gibbons et al. 
(Gi59). 


4. Neutron escape from MnSO, bath. 


4 
of the NBS standard source. This correction was computed by dividing 


Neutron escape from the MnSO, bath was significant only in the case 
the outer boundary of the bath into zones, computing the neutron flux in 
each surface zone, and weighting by the solid angle subtended at the center 
of the bath. The fraction f of neutrons escaping in each zone was com- 
puted by the formula 


f = 2 2°]. soumeen (Ge59) 
? 


where r is the distance from the source to the edge of the tank in 


centimeters. 


d/r 
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